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Abstract 

By performing fully 3D symmetry-unrestricted Skyrme-Hartree-Fock- 
Bogoliubov calculations, we discuss shape coexistence and possibility of 
exotic deformations simultaneously breaking the reflection and axial 
symmetries in proton-rich N = Z nuclei: 64 Ge, 68 Se, 72 Kr, 76 Sr, 80 Zr 
and 84 Mo. Results of calculation indicate that the oblate ground state of 
68 Se is extremely soft against the I33 triangular deformation, and that 
the low-lying spherical minimum coexisting with the prolate ground 
state in 80 Zr is extremely soft against the I32 tetrahedral deformation. 
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1 Introduction 

The Hartree-Fock-Bogoliubov (HFB) method with the Skyrme interactions is 
one of the standard approaches in nuclear structure research [ILQ. In the last 
two decades it has become possible to solve the HFB equations directly in the 
coordinate mesh space ||f§]. In recent years, in order to investigate the struc- 
ture of drip-line nuclei, the need for such coordinate-space HFB calculations 
has been greatly increased and intensive analyses have been made for neutron 
radii and skins in spherical neutron- rich nuclei |5|-11|: Since the easier HF 



plus BCS method breaks down when treating the pairing correlation in weakly 
bound systems due to a leakage of nucleons into the continuum, we need to 
calculate the mean-field (particle-hole) correlations and the pairing (particle- 
particle) correlations selfconsistently in order to preserve confinement of the 
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nuclear density while allowing the pairing excitations to positive energy reso- 
nant states || (see, e.g. |12] and references therein for mean-field approaches 
other than the Skyrme-HFB method). 

Recently, Terasaki, Heenen, Flocard and Bonche [fL3, 14| have removed the 
restriction of spherical symmetry in solving the coordinate-space Skyrme-HFB 
equations in order to investigate the possibility to get three-dimensional (3D) 
deformed solutions in neutron rich nuclei. In their works, a Skyrme interaction 
is used to describe the Hartree-Fock (HF) Hamiltonian while a density depen- 
dent zero-range interaction is used for the pairing channel. The mean-field 
HF equations are solved by the imaginary-time evolution method [15|] in a 3D 
cubic mesh space while the HFB equations are solved in terms of the two-basis 
method developed earlier in IIB,0. The discretization in 3D mesh space has 
the advantage over methods relying on an expansion in the harmonic-oscillator 
basis that nuclei with exotic deformations can be treated at the same level of 
accuracy fll~8|-f2~0|j . In these works, however, reflection symmetries with respect 
to three planes are imposed for the nuclear density so that only one spatial 
octant is needed to solve the HFB equations. 

The major purpose of this paper is to extend their method by removing 
the symmetry restrictions mentioned above and investigate the possibility of 
exotic shapes simultaneously breaking the axial and reflection symmetries in 
the mean field. For this purpose, we have constructed a new computer code 
that carries out Skyrme-HFB calculations in the 3D Cartesian-mesh space 
without imposing any restrictions on the spatial symmetry. Recently, on the 
basis of the Skyrme HF plus BCS calculations with no restriction on the nuclear 
shape, Takami, Yabana and Matsuo [|2T],[£J suggested that the oblate ground 
state of 68 Se is extremely soft against the Y33 triangular deformation, and 
that the low-lying "spherical" minimum coexisting with the prolate ground 
state in 80 Zr has the Y32 tetrahedral shape. As the first application of a fully 
3D, symmetry- unrestricted Skyrme HFB method with the use of the density- 

, H|, |TT|, P3"|-|3Ti| , we investigate in 



dependent, zero-range pairing interaction 
this paper shape coexistence and possibility of non-axial octupole deformations 
in proton rich N = Z nuclei in the A = 64 — 84 region and examine the above 
predictions. These nuclei are especially interesting objects to study, since 
proton and neutron deformed shell effects act coherently and rich possibilities 
arise for coexistence and competition of different shapes (see [[H] for earlier 
references). In recent years, active experimental studies of these nuclei are 
going on by means of combinations of radioactive nuclear beams and new 
gamma-ray and charged particle detector systems (see |]32|-p5|1 for reviews). It 



should be noted here that, although extensive theoretical calculations and rich 
experimental evidences have been accumulated for axially symmetric octupole 
(F30) deformations, as reviewed in |35|,[J7|], only a few calculations using Woods- 
Saxon Strutinsky methods are available p8|-Bl|] except for light nuclei, and no 



firm experimental evidence exists up to now concerning the non-axial octupole 
(Y31, Y32, Y33) deformations in the mean fields. For light nuclei, non-axial 
octupole deformations have been discussed [|2|-446| in connection with alpha- 
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cluster structures P7| ; for instance, a triangular structure of 12 C p£2|,[4ll| and 
a tetrahedral shape for 16 [^5|,|46|. 



Our motive for developing the coordinate-space Skyrme-HFB method is 
not only to investigate the possibilities of emergence of new types of symmetry 
breakdown in the ground states of proton-rich and neutron-rich nuclei, but 
also to investigate, in the future, low-lying modes of excitation of such unsta- 
ble nuclei by means of the RPA and the Selfconsistent Collective Coordinate 
(SCC) method [48| on the basis of the HFB basis thus obtained. We intend to 
proceed in parallel with other calculations with the use of more phenomeno- 
logical shell model potentials and separable interactions. The Skryme-HFB 
method is suited for this aim, as it provides a local mean-field potential so 
that such a comparative study is easy. 

In Section 2, a brief account of the method of the coordinate-space Skyrme- 
HFB calculation is given. In Section 3, results of numerical calculation are 
presented and discussed. In Section 4, a conclusion is given. 



2 Skyrme-HFB calculation 

2.1 Two basis method 



For convenience, we here recapitulate the two basis method fl3|, |14|, |16|, [17 



adopted as the algorithm of our computer code. In this method, the imaginary- 
time evolution method is combined with a diagonalization of the HFB Hamil- 
tonian matrix to construct the canonical basis. 

We first determine the single-particle wave functions <pi satisfying the HF 
equations 

h\p(r)]<l> i {r)=e i Mr) (1) 



by means of the imaginary-time evolution method [TjJ. Here h, £, and p(r) 



denote the mean-field Hamiltonian, the single-particle energies and the total 
nuclear density, respectively. (The isospin index r is omitted for simplicity.) 
We next diagonalize the HFB Hamiltonian matrix [0 



h-X A \ fU k \ „ [U k 



A * - h * + X J\V k J- Ek [v k ) (2) 
to get the one-body density matrix p and the pairing tensor k: 



p = V*V T , k = V*U T . (3) 

We then diagonalize the density matrix p and obtain the occupation coefficients 
n a and the unitary transformation W which relates the HF wave functions 4>i 
to the canonical basis wave functions ip a : 
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PH = 'En a W ka Wl (4) 



^«=E^i«' ( 5 ) 

j 

In the canonical basis ip a , the HFB density matrix in the coordinate space is 
diagonal: 

p(r,r') =J2 n afa(r)(p a (r')*. (6) 

a 

These steps are repeated until the convergence is achieved. 

The single-particle wave functions and densities are represented on a full 3D 
Cartesian mesh space within a spherical container. In the present calculation, 
the radius of the spherical container and mesh spacing are set to R mes h = 
10.0 fm and h = 1.0 fm, respectively. Tajima et al. have carefully 

examined possible errors due to the use of the mesh size h = 1.0 fm and they 
found that, since discretization errors are essentially independent of the nuclear 
shape, deformation energies obtained with this mesh size are quite accurate 
(see also 0]). Actually, we have constructed the new Skyrme-HFB code by 



extending the cranked Skyrme-HF code |52[ written previously and applied to 
the investigation of the yrast structure of 32 S, so that the cranking term can 
be included. In this paper, however, we examine only the cases of zero angular 
momentum. 

2.2 The Skyrme plus density- dependent pairing interactions 

We use the SIII parameter set fl53j of the Skyrme interaction for the mean-field 
(particle-hole) channel, which has been successful in describing systematically 
the ground state quadrupole deformations in proton and neutron rich Kr, Sr, 
Zr and Mo isotopes [jl9[ and in a wide area of nuclear chart [|{|. For the pairing 



(particle-particle) channel, we use the density- dependent zero-range interaction 
[|i~3l, |l4| [l7|. |23"H30l , which has been successful in describing, for instance, the 
odd-even staggering effects in charge radii, 



V pair (n, r 2 ) = y (l - P a ) (l - ^\ 5 (n - r a ) (7) 



with the notation of ||17|| , where the strength Vq and the density p c are param- 
eters and P a denotes the spin exchange operator. For these parameters, we 
use the standard values |4],[17|: V Q = -1000.0 MeV-fm 3 , p c = 0.16 fm" 3 . The 
pairing interaction is smoothly cut off at 5 MeV above the Fermi energy in the 
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same way as in [ I9[ . For a more general form of the density dependent pairing 
interaction, we refer pi| . [55| . 

To check the dependence on the Skyrme-interaction parameter sets, we 
make calculations with the SkM* |56| and SLy4 |57j sets for an example of 68 Se. 



We refer to a recent work by Reinhard et al. ]5lJ for a detailed and systematic 
study of shape coexistence phenomena in relation to the properties of various 
versions of the Skyrme interaction. We shall also check the dependence on the 
pairing strength V adopted. 

2.3 Constrained HFB calculation 

In order to investigate the deformation properties away from the HFB equi- 
librium points, we perform constrained HFB calculations with the use of the 
quadratic constraints for the mass-quadrupole (octupole) moments |59| to ob- 




tain the energy surfaces as functions of the quadrupole (octupole) deforma- 
tions. Because no spatial symmetry is imposed on the 3D mesh space, the 
center of mass and the directions of the principal axes of the nucleus can move 
freely without affecting the total energy. To evaluate the physical quantities 
like deformation parameters, it is crucially important to fulfill the constraints 
to keep the center of mass, 



and the directions of the principal axes, 



(Eta/)*) = (£M<) = (EM*) = o. (9) 

These requirements are taken care of by means of the quadrupole con- 
straints on these conditions as in our previous study JEJ. 

2.4 Deformation parameters 

As measures of the deformation, we calculate the mass-multipole moments, 



An f 

a lm = J r l X lm (tt) p(r)dr, (m = -/,••■,/) (10) 

where p{r) = £^ \<p a (r)| 2 , R = 1.2A 1 / 3 fm and X lm are real bases of the 

a 

spherical harmonics, 
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Xio = Y l0 , (11) 
1 

V2' 



X l\m\ = —fc(Yl-\m\+Yl-\ m \), (12) 



Xl-\m\ = -|(1/ H -^h)- ( 13 ) 

Here the quantization axis is chosen as the largest (smallest) principal axis for 
prolate (oblate) solutions. We then define the quadrupole deformation param- 
eter /?2, the triaxial deformation parameter 7, and the octupole deformation 
parameters /?3 and (3^ m by 



"20 = P2 cos 7, a 22 = P2 sin 7, (14) 
/ 3 \ 1/2 

& = E «L , ft™ = («L + «s- m ) V2 (m = 0, 1, 2, 3) . (15) 

\m=-3 J 

For convenience, we also use the familiar notation — fa for oblate shapes 
with (#2,7 = 60°). 



3 Results and discussions 

3.1 Quadrupole deformations 

The solutions of the Skyrme-HFB equations obtained in the numerical calcu- 
lations for 64 Ge, 68 Se, 72 Kr, 76 Sr, 80 Zr and Mo are summarized in Table 1. 
The calculated ground-state shape changes from triaxial ( 64 Ge), oblate ( 68 Se, 
72 Kr), large prolate ( 76 Sr, 80 Zr), to spherical shape ( 84 Mo) with increasing 
N(= Z). For 68 Se, 72 Kr, 76 Sr, 80 Zr and 84 Mo, we obtain two or three local 
minima close in energy, indicating shape coexistence. These gross features are 



consistent with available experimental data [BC-SB| and previous theoretical 
calculations 0,^0,(11,0^. 

The potential energy curves obtained by the constrained HFB calculations 
are displayed in Fig. [I] as functions of the quadrupole deformation parameter 
fa and in Fig. |] as functions of the triaxial deformation parameter 7. Below 
we remark on some specific points. 

As seen in Fig. [|, the calculated potential energy curve for 64 Ge is rather 
shallow with respect to the 7 degree of freedom so that this nucleus may 
be regarded as "7-soft." This result is consistent with the experimental in- 



dication |K§ and also with the shell model calculation by the Monte Carlo 



diagonalization method [75]. 

Quite recently, an excited prolate band coexisting with the ground-state 
oblate band has been found in 68 Se Their quadrupole deformations are 
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estimated as fa ~ 0.27 and fa ~ —0.27, respectively. Although the prolate 
excited band-head + state has not yet been observed, its excitation energy is 
estimated to be about 0.6 MeV. Our calculated energy difference between the 
prolate and the oblate HFB solutions, 0.52 MeV, is in good agreement with 
this experimental data. The barrier between the prolate and the oblate minima 
is about 3 MeV in the plot with respect to fa in Fig. [I], but it is only about 0.3 
MeV in the plot with respect to the triaxial deformation parameter 7 in Fig|| 
It might be considered that, if the barrier is so low, the two bands built on the 
prolate and the oblate solutions interact strongly so that the shape coexistence 
picture is too much perturbed in contradiction with the experiment In 



our view, however, description of dynamics by going beyond the static mean- 
field approximation is necessary in order to discuss the interaction between 
the oblate and the prolate structures. In any case, understanding this shape 
coexistence dynamics is an interesting subject for future. 

The second minimum with fa »s 0.66 seen in the potential energy curve for 
84 Mo in Fig. |T| may be regarded as a superdeformed solution, since it is related 
to the Z = N = 42 deformed shell gap |68| formed by occupying the down- 
sloping [431] 1/2 levels from the upper major shell by two protons and two 
neutrons. This second minimum was also obtained in f2~I |. It offers an inter- 



esting possibility that a superdeformed rotational band might be observed at 
such a low excitation energy as about 1.5 MeV. From a viewpoint of deformed 
shell structure, the ground-state solutions for 76 Sr and 80 Zr have character- 
istics different from the second minimum in 84 Mo and may be distinguished 
from the superdeformation, although they have large prolate deformations of 
fa « 0.5. 

3.2 Non-axial octupole deformations 

As a result of the Skyrme-HFB calculations for proton-rich N = Z nuclei 
from 64 Ge to 84 Mo (summarized in Table 1), we have found equilibrium shapes 
with finite non-axial octupole deformations for 68 Se and 80 Zr. The density 
distribution at the HFB local minimum for 68 Se with the triangular deforma- 
tion superposed on the oblate shape and that for 80 Zr with the tetrahedral 
deformation are illustrated in Fig. |3|. 



In addition to the two cases mentioned above, Takami et al. 211 and Matsuo 



et al. |22] obtained, in their Skyrme-HF plus BCS calculations, finite equilib- 
rium values of octupole deformations superposed on an oblate shape in 76 Sr 
and also on a near spherical shape in 84 Mo. According to their calculations, 
the potential-energy curves are very soft with respect to the octupole deforma- 
tion degrees of freedom especially in the four cases mentioned above. In order 
to see the properties of the potential-energy curve in the neighborhood of the 
HFB equilibrium points and to make a better comparison with the results of 
Refs. J2~T|, [Z2"| , we have carried out constrained HFB calculations with respect 
to the fa m (m = 0, 1, 2, 3) degrees of freedom about the local minima (seen in 
Fig. 1) of the quadrupole deformation energies. 
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Figures 4, 5 and 6 show the potential-energy curves with respect to the 
octupole deformation parameters /3 3m about the oblate, the spherical and the 
prolate (or triaxial) minima of the quadrupole deformation energy curves, re- 
spectively. These curves are obtained by the constrained HFB calculations 
with the octupole operators r 3 X 3 | m | as constraints. We see that the oblate 
shape of 68 Se is extremely soft against the triangular (^33 ) deformation and 
that the spherical shape of 80 Zr is extremely soft against the tetrahedral (P32) 
deformation, in agreement with those of the Skyrme-HF plus BCS calculations 
of Refs. [21, 22] . The oblate shape of 76 Sr is fairly soft with respect to the ^32 



and /?33 deformations and the spherical ground state of 84 Mo is barely stable 
against all /3 3m degrees of freedom, especially against /3 30 . In [|^] an oblate so- 
lution with a finite equilibrium value of ^32 is obtained for 76 Sr, while a similar 
solution for 76 Sr but with a finite equilibrium value of ^33 and also a nearly 
spherical solution for 84 Mo with a finite equilibrium value of (3%$ is reported 
f2~T| . Although such details differ depending on the treatment of the pairing 
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correlations, the basic features, i.e., the softness to both /5 32 and ^33 of the 
oblate shape of 76 Sr and the softess to /5 30 of the spherical shape of 84 Mo are 
in common between the present HFB calculations and those of ||2l],22]. Gener- 
ally speaking, Figs. 4-6 indicate that the oblate shapes are softer for octupole 
deformations /3 3m with higher values of m, while the prolate shapes favor lower 
values of m. 

Below we focus our attention on the triangular deformation in 68 Se and the 
tetrahedral shape in 80 Zr and discuss about the microscopic origins of them. 



Triangular deformation in 68 Se 



Generally speaking, octupole correlations are associated with strong cou- 
plings between the shell- model orbits with Al = Aj = 3 |36| , |37f . In the 
A = 64 — 84 region under consideration, they are lgg/2 and 2p 3 / 2 . In order 
to understand why the oblate shape in 68 Se is unstable (or extremely soft) 
against the triangular deformation, however, we need to examine the interplay 
of the quadrupole and octupole deformation effects. Namely, as explained be- 
low, the emergence of the triangular deformation is strongly correlated with 
the magnitude of the oblate deformation. 

When 68 Se (N = Z = 34) is oblately deformed, the high tt levels [404] § 
and [413] I stemming from the lgg/2 orbit go down in energy and approach 
the Fermi surfaces for N = Z = 34 and strong Y" 3 3 couplings with [301] | and 
[310] I levels (associated with the 2p 3 / 2 orbit) take place. These Y33 coupling 
effects are seen as repulsions between these levels in Fig. [7| which displays 
the neutron single-particle energies as functions of the triangular deformation 
parameter j3^. Here, the single-particle energies mean eigenvalues of the HF 
Hamiltonian with the density p (r) determined by the HFB equations, and the 
asymptotic Nilsson quantum numbers are used only for convenience of labeling 
these levels: they are, of course, not good quantum numbers. 

In this figure, results of calculation with use of the SkM* and SLy4 interac- 
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tions are also shown for comparison. We note that the Y33 coupling effects are 
slightly weaker in the case of the SkM* and SLy4 interactions in comparison 
with the case of the SIII interaction. This is because the spacings between 
the levels coupled by the Y33 operator are the smallest for the SIII interaction: 
The spacings at the oblate equilibrium deformations between the [404] 9/2 and 
[301]3/2 levels are about 2.8, 3.4 and 3.6 MeV, and those between the [413]7/2 
and [310] 1/2 levels are about 3.8, 4.1 and 4.2 MeV for the SIII, SkM* and 
SLy4 interactions, respectively. Thus, as shown in Fig. ^, the potential energy 
curve with respect to the triangular /3 33 deformation is softest for the case 
of the SIII interaction, although they are soft also for the cases of the SkM* 
and SLy4 interactions. Note that, in making this comparison, we have chosen 
the pairing-interaction strength Vq such that the resulting pairing gaps A take 
about the same values for calculations with different Skyrme interactions (in 
order to make the effects of the pairing correlations approximately the same 
for all cases), as shown in the right-hand part of Fig. [8[ 

The importance of the triangular Y33 deformation superposed on the oblate 
shape was previously pointed out by Frisk, Hamamoto and May JFHJ in terms of 
a two-level model as well as the modified oscillator model which simulates the 
one-particle spectra in an infinite-well potential. Our result of the Skyrme- 
HFB calculation provides a realistic example which is consistent with their 
arguments. 

Tetrahedral deformation in 80 Zr 

As shown by Hamamoto, Mottelson, Xie and Zhang [|79|| , the tetrahedral 
symmetry associated with the Y 32 deformation brings about a bunching of 
the single-particle levels and create a remarkable shell structure: The iV = 
Z = 40 is one of the magic numbers for such tetrahedral shapes. Such a shell 
effect is common to various finite Fermion systems, and in fact the tetrahedral 
deformation has been predicted, for instance, for sodium clusters consisting of 
40 atoms by the density functional Kohn-Sham calculation 0, [Hp, in which 
there is no spin-orbit coupling. The instability of the spherical shape of 80 Zr 
against the Y32 deformation, as exhibited in Fig. [5|, is evidently connected to 
the magic number N = Z = 40 for the tetrahedral shape. 

Figure ^| shows the single-particle energy diagrams as function of octupole 
deformation parameter /3 3rn (m = 0, 1, 2, 3). As expected, we can see for the case 
of m = 2 a remarkable bunching of single-particle levels and an increase of the 
shell gap at N = 40 with increasing ^32, while the other octupole deformations 
(m = 0, 1, 3) do not exhibit such a feature. Looking into details, one notices a 
fine splitting of the lgg/2 level into three levels which correspond to irreducible 
representations of the double tetrahedral (spinor-T^) group M,H5[|; a twofold- 
degenerate level and two fourfold-degenerate levels. 

Thus, the tetrahedral shell gap at N = Z = 40 emerges even under the 
presence of the strong spin-orbit coupling. It should also be noted that the 
tetrahedral minimum is obtained in the calculation selfconsistently including 
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the pairing correlations. 



3.3 Pairing gaps 

In this subsection, we first examine dependence of the pairing gaps on defor- 
mations, and then discuss dependence of the non-axial octupole deformations 
on the pairing strength. The result of calculation for the pairing gaps at equi- 
librium deformations in each nucleus is listed in Table 1. As the pairing gaps 
in the HFB theory depend on single-particle levels, the numbers listed in this 
Table are averages of the diagonal elements in the HF basis, over 5 MeV 
interval in the vicinity of the Fermi surfaces. 

In the literatures, slightly different quantities like averages of the diagonal 
matrix elements in the canonical basis, A adl , weighted by the coefficients of 
the Bogoliubov transformation, u a v a l82|-|34"f or ^ ||, are used for similar 



purposes. Figure |TU] compares these quantities for the case of triangular defor- 
mations superposed on the oblate shape in 68 Se. We see that the two average 
quantities, (A-) and (A a aU a v a ) , are approximatelty equal. We also confirm 
that the averages do not significantly depend on the averaging interval. 



Figures [11], [12| and [13] display the variation of the pairing gaps with the 
quadrupole deformation parameter fa, the triaxial deformation parameter 7, 
and the octupole deformation parameters (3z m {m = 0,1,2,3), respectively. 
We observe that gross features of deformation dependence of the pairing gap 
correlate with the corresponding potential-energy curves displayed in Figs. 
1, 2 and 4-6. Such correlations are rather easy to be understood from the 
behavior of the single-particle level density near the Fermi surface, i.e., from 
the well-known (spherical or deformed) shell effects that the level density near 
the Fermi surface becomes relatively low in the vicinities of the local minima of 
the potential energy curve |^J. Thus, the pairing correlation becomes weaker 
and the paring gap decreases near the local minima. On the other hand, the 
level density becomes relatively high and the pairing gap increases near the 
local maxima of the potential-energy curve. 

Because of significant shape changes in the sequence of isotopes (isotones) 
in the A = 64 — 84 region, it is not always easy to extract the magnitudes 
of pairing correlations from experimental odd-even staggerings of binding en- 
ergies and to assess the appropriateness of the pairing-interaction strength 
Vq = —1000 MeV- fm 3 used in our HFB calculations. Quite recently, how- 



ever, Satula, Dobaczewski and Nazarewicz |86| have proposed a method for 
separating out the pairing correlation effects from the deformed mean-field 
(single-particle energy) effects on the odd-even staggerings, and evaluated av- 
erage pairing gaps; these are in the range 1.0 — 1.6 MeV for the mass region 
under consideration [37]] . We note that these values agree rather well with the 



well-known global trend A = 12/ v A MeV f88| , which are in the range 1.3 — 1.5 
MeV for A = 64 — 84. Our calculated values of the pairing gaps, listed in Table 
1 and drawn in Figs. |TTLff2lff3L mostly lie in this range of values, so that we 



may say that the adopted strength for V is reasonable. 
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Another possible source of ambiguity in evaluating the pairing gaps is the 
proton-neutron isoscalar pairings which are expected to play an important 



role in the N = Z nuclei (see, for example, [89,9(| and references therein). 
We have assumed that such isoscalar pairings are absent in the states under 
consideration. Although this assumption should be examined, there are some 



experimental indications |p3|,|91|j that this may be a fairly good approxima- 
tion. It is clear that we need a more systematic and detailed investigations, 
both theoretical calculations and experimental explorations, for a better un- 
derstanding of the pairing correlations in the proton-rich N = Z nuclei in the 
A = 64 — 84 region. 

In order to examine the sensitivity of the calculated results to the strength 
Vq of the pairing interaction, we have made a calculation of the potential energy 
curve about the oblate shape in 68 Se as a function of the triangular octupole 
deformation parameter /3 33 for Vq = —900,-1000 and —1100 MeV-fm 3 . The 
result is shown in Fig. [14]. As expected, the potential energy curve becomes 
shallower with increasing (absolute value of) Vq. Thus, the local minimum at 
/3 33 ~ 0.10 disappears with 10% increase of the (absolute) value of Vq. In any 
case, the potential is so shallow that we cannot associate a definite physical 
significance with the equilibrium values of /3 33 . We can still draw from these 
calculations an important conclusion that the oblate ground state of 68 Se is 
extremely soft with respect to the triangular octupole deformation. 

3.4 Discussions 

Actually, we need a more detailed investigation on the physical implication of 
the extremely soft potentials like those with respect to the triangular deforma- 
tion in 68 Se and for the tetrahedral shape degree of freedom in 80 Zr. As is well 



known in the case of the axially symmetric Y 30 octupole deformation [92-96], 
a definite minimum develops at finite value of /3 3 o after the parity projection 
when the mean-field potential is very soft with respect to /3 3 o. For the case of 
non-axial octupole deformations, a similar effect of the parity projection has 



been demonstrated by Takami, Yabana and Ikeda [|2| for light nuclei. It re- 
mains to be examined whether or not the situation is similar for the non-axial 
octupole deformations in medium-mass nuclei under consideration. 

More generally speaking, investigations of modes of excitation and of exci- 
tation spectra associated with the instabilities toward the non-axial octupole 
shape deformations is one of the major challenges for future. The present 
paper should be regarded as providing a HFB mean-field basis for a study 
of dynamics by means of methods like the quasiparticle RPA and the SCC 
method 
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4 Conclusion 



We have constructed a new computer code that carries out Skyrme-HFB cal- 
culations in the 3-dimensional Cartesian-mesh space without imposing any 
restriction on the spatial symmetry, and investigated shape coexistence and 
non-axial octupole deformations in proton-rich N = Z nuclei, 64 Ge, 68 Se, 72 Kr, 
76 Sr, 80 Zr and 84 Mo. The ground state shape changes from triaxial ( 64 Ge), 
oblate ( 68 Se, 72 Kr), large prolate ( 76 Sr, 80 Zr), to spherical ( 84 Mo) as N(= Z) 
increases, in agreement with the available experimental data and the previous 
theoretical calculations. The extreme softness toward the Y33 triangular de- 
formation of the oblate ground state of 68 Se and that toward Y 32 tetrahedral 
deformation of the excited spherical minimum of 80 Zr, pointed out by Takami 
et al. pi| , |2"2]| on the basis of the Skyrme-HF plus BCS calculations, have been 
confirmed by the fully selfconsistent Skyrme-HFB calculations with the use of 
the density-dependent zero-range pairing interaction. 

The symmetry-unrestricted Skyrme-HFB computer code constructed in 
this work provides a selfconsistent mean-field basis for future investigation 
of collective modes of excitation in neutron-rich nuclei with neutron skins as 
well as in proton-rich nuclei. 
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Oblate Spherical Prolate 


64 Ge 


g.s. 

/3,7 = 0.27,25° (triaxial) 
A» = 0.0 
A p = 1.25, A n = 1.12 


68 Se 


g.s. 0.52 
7 = 0.28,60° 0, 7 = 0.26,0° 
& = Paa « 0.08 & = 0.0 
A p = 1.28, A n = 1.13 A p = 1.29, A n = 1.15 


72 Kr 


g.s. 0.92 
/3, 7 = 0.32,60° /3, 7 = 0.40,0° 
P 3 = 0.0 /3 3 = 0.0 
Ap = 1.03, A n = 1.23 Ap = 1.25, A n = 0.92 


76 Sr 


1.79 g.s. 
P, 7 = 0.30,60° /3, 7 =0.51,0° 
Ps = P33 « 0.0 ft = 0.0 
Ap = 1.47, A n = 1.43 Ap = 0.67, A n = 0.50 


80 Zr 


0.86 1.01 g.s. 
P, 7 = 0.20,60° P, 7 = 0.0,0° P, 7 = 0.51,0° 
& = 0.0 p 3 = p 32 w 0.15 /5s = 0.0 
Ap = 1.02, A n = 0.82 A p = 0.68, A n = 0.39 A p = 0.79, A„ = 0.78 


84 Mo 


0.20 g.s. 1.52 
P, 7 = 0.16,60° 7 = 0.0,0° /3, 7 = 0.66,0° 

A, = 0.0 p 3 = p 30 w 0.0 & = 0.0 
Ap = 1.46, A n = 1.42 Ap = 0.74, A n = 0.72 A p = 0.0, A„ = 0.0 



Table 1: Solutions of the HFB equations for proton-rich N = Z nuclei in the 
A = 64 — 84 region. For each nucleus, numbers in the first line indicate excitation 
energies measured from the ground state. The symbol ~ indicates that the potential 
energy curve is extremely shallow about the equilibrium value. Pairing gaps A p and 
A n are here defined as averages of diagonal elements over 5 MeV interval around 
the Fermi surface, and their values (in MeV) at the equilibrium deformations are 
listed. 
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Figure 1: Potential energy curves calculated by the constrained Skyrme-HFB 
method for 64 Ge, 68 Se, 72 Kr, 76 Sr, 80 Zr and 84 Mo are drawn as functions of the 
quadrupole deformation parameter /?2. The SIII interaction is used for the particle- 
hole channel, while the density-dependent pairing interaction with Vq = —1000.0 
MeV-fm 3 and p c = 0.16 fm~ 3 is used for the particle-particle channel. 
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Figure 2: Potential energy curves calculated at fixed by the constrained Skyrme- 
HFB method for 64 Ge, 68 Se, 72 Kr, 76 Sr, 80 Zr and 84 Mo are drawn as functions of 
the triaxial deformation parameter 7. The effective interactions used are the same 
as in Fig. 
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Figure 3: Density contour surfaces at the half central density of the Skyrme- 
HFB solution with the oblate plus triangular shape (/% = —0.28, /?33 = 0.08) for 
68 Se (left-hand side) and that with the tetrahedral shape (fo = 0.00, P32 = 0.15) for 
80 Zr (right-hand side), listed in Table 1. 
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Figure 4: Potential energy curves calculated by the constrained Skyrme-HFB 
method are drawn as functions of the octupole deformation parameters /33 m (m = 
0,1,2,3) about the oblate minima (seen in Fig. 1) of the quadrupole deformation 
energies. One of the f3^ m {m = 0, 1, 2, 3) is varied while the other /?3 m 's are fixed to 
zero. The effective interactions used are the same as in Fig. [l]. 
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Figure 5: Same as Fig. |] but about the spherical minima. 
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ure 6: Same as Fig. |] but about the prolate minima (the triaxial minimum in 
case of 64 Ge). 
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-16 t A 

0.05 0.1 0.15 0.2 0.25 



Figure 7: Neutron single-particle energies for 68 Se plotted as functions of the oc- 
tupole deformation parameter (3^ about the oblate shape. Here, the single-particle 
energies mean eigenvalues of the HF Hamiltonian with the density p (r) determined 
by the HFB equations. Results for the SIII, SkM* and SLy4 parameter sets are com- 
pared. Equilibrium quadrupole deformations obtained for each Skyrme interaction 
are (3 2 = -0.28, -0.25 and -0.24 for SIII, SkM* and SLy4, respectively. Solid (bro- 
ken) lines indicate levels which have positive (negative) parity in the limit ^33 = 0. 
The projection of the angular momentum on the symmetry axis, f2, is a good quan- 
tum number only at ^33 = 0. The arrows indicate the AO, = 3 coupling associated 
with the triangular Y33 deformation as discussed in the text. The single-particle 
spectrum for protons is almost the same as for neutrons. 
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Figure 8: Comparison of the HFB potential energy curves for 68 Se about the 
oblate shape as functions of the triangular deformation parameter (5^ , calculated for 
different versions of the Skyrme interaction (left-hand side) . The pairing-interaction 
strengths Vq are chosen such that the average pairing gaps become approximately 
equal for all Skyrme interactions (as displayed in the right-hand side) . The calculated 
deformation parameter /3 2 are -0.28, -0.25 and -0.24 for the SIII, SkM* and SLy4 
interactions, respectively. 
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Figure 9: Neutron single-particle energies for 80 Zr plotted as functions of the 
octupole deformation parameters (3z m (m = 0,1,2,3) about the spherical shape. 
Here, the single-particle energies mean eigenvalues of the HF Hamiltonian with 
the density p (r) determined by the HFB equations. The SIII interaction is used. 
Solid (broken) lines indicate levels which have positive (negative) parity in the limit 
/?32 = 0. The single-particle spectrum for protons is almost the same as for neutrons. 
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Figure 10: Comparison of differently defined average pairing gaps for 68 Se, 
plotted as functions of the triangular deformation parameter ^33 superposed on 
the oblate shape. Here, (A U ) AE = J2i fifi9iSrAii/T,i fifi9i9i, (uvA) cano = 
{(f a |A| (f a ) /J2 and (v 2 A) cano = J2 a «S (<Pa I A l <Pa) I E« v 2 a 

@, where /< = (1 + exp[( £i - \ F - AE/2)/ fj])' 1 ^, 9i = (1 + exp[( £i - X F + 
A£;/2)//i])- 1 / 4 with AE = 3 or 5 MeV. 
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Figure 11: Variations of the pairing gaps A T (r= p, n) calculated by the con- 
strained Skyrme-HFB method as functions of the quadrupole deformation parame- 
ter /?2 for 64 Ge, 68 Se, 72 Kr, 76 Sr, 80 Zr and 84 Mo. The effective interactions used are 
the same as in Fig. ||. 
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Figure 12: Variations of the pairing gaps A r (r = p, n) calculated by the con- 
strained Skyrme-HFB method as functions of the triaxial deformation parameter 7 
at fixed P2 for 64 Ge, 68 Se, 72 Kr, 76 Sr and 80 Zr. The effective interactions used are 
the same as in Fig. [l|. 
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Figure 13: Variations of the pairing gaps A T (r= p, n) calculated by the con- 
strained Skyrme-HFB method as functions of the octupole deformation parameter 
P3 m (m = 0, 1, 2, 3) about the local minima (seen in Fig. 1) of the quadrupole defor- 
mation energies for 64 Ge, 68 Se, Kr, 76 Sr, 80 Zr and 84 Mo. The effective interactions 
used are the same as in Fig. [jj 
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Figure 14: Comparison of the potential energy curves (left-hand side) and average 
pairing gaps for protons (right-hand side) calculated by the constrained Skyrme- 
HFB method as functions of the triangular deformation parameter ^33 about the 
oblate shape for 68 Se with use of different strengths Vq of the density-dependent 
pairing interaction (and with the same SIII interaction ). 
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